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Abstract 



We perform a calculation of Sivers function in a light-cone SU(6) quark-diquark 
t"*- ■ model with both scalar diquark and vector diquark spectators. We derive the trans- 

, verse momentum dependent light-cone wave function of the proton by taking into 

f— ^ I account the Melosh-Wigner rotation. By adopting one-gluon exchange, we obtain 

■^ ' a non-vanishing Sivers function of down quark from interference of proton spin 

^^ amplitudes. We analyze the —tt^ weighted Sivers asymmetries in 7r+, vr^ and vr^ 



i-i^ I electroproduction off transverse polarized proton target, averaged and not averaged 

by the kinematics of HERMES experiment. 
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1 Introduction 



The investigation of single-spin asymmetries in hadronic reaction has attracted 
great attention because it can uncover some novel structure of the nucleon. 
Various theoretical explanations have been put forward. One of them is the so- 
called Sivers effect [1], which indicates that the transverse motion of quarks in 
transverse polarized proton can lead to an azimuthal production asymmetry. 
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The effect can be factorized as a convolution of a T-odd transverse momen- 
tum dependent parton distribution function (Sivers function) /f!^(x,k^), the 
number density of unpolarized quarks in transversely polarized proton, with 
the unpolarized fragmentation function. However, for several years the Sivers 
function was believed to vanish due to the time reversal invariance property of 
QCD. An alternative mechanism hereafter named the Collins effect was pro- 
posed in Ref. [2]: the chiral-odd fragmentation function convoluted with the 
transversity distributions of the proton can generate the single-spin production 
asymmetries. Collins' idea motivated the experimental and theoretical efforts 
to reveal the transversity properties of the nucleon [3]. Semi-inclusive DIS 
(SIDIS) experiments on longitudinally polarized target by HERMES collabo- 
ration [4] and on transversely polarized target by SMC collaboration [5] showed 
significant azimuthal asymmetries. During the same time, a number of theo- 
retical analyses on these observed asymmetries have been undergone in terms 
of the Collins effect [6,7,8,9]. Recently, Brodsky, Hwang and Schmidt [10,11] 
proposed a new mechanism of producing single spin asymmetries, and they 
demonstrated that the final state interactions from gluon exchange between 
the outgoing quark and target spectator system can lead to single-spin asym- 
metries in SIDIS or Dell-Yan process at leading-twist level in perturbative 
QCD. It was shown soon by Collins [12] that this mechanism can be recog- 
nized as Sivers effect, and the Sivers function is not vanish due to the presence 
of Wilson line in the operator definitions of parton density. This implies that 
the gauge link term in the gauge-invariant definition of the parton distribu- 
tion function can generate the final or initial state interactions [13], which are 
necessary to produce a phase difference for a nonzero Sivers function. 

The allowance of T-odd parton distribution opens up a wide range of phe- 
nomenological applications. In this work we attempt to give a full analysis 
of Sivers asymmetries in SIDIS in the scope of final state interactions. Based 
on the quark-scalar-diquark model employed in [10], we incorporate vector 
diquark structure into the proton wave function through a simple relativistic 
quark spectator-diquark model which is formulated in the light-cone frame. 
This model was originally proposed in order to study deep inelastic lepton 
nucleon scattering [14,15], based on the quark-parton model picture that deep 
inelastic scattering is well described by the impulse approximation, in which 
the incident lepton scatters incoherently off a quark in the nucleon, with the 
remaining nucleon constituents treated as a quasiparticle spectator. After tak- 
ing into account Melosh rotation effects, this model is in good agreement with 
the experimental data of polarized deep inelastic scattering, and the mass dif- 
ference between the scalar and vector spectators reproduces the up and down 
valence quark asymmetry [16]. By advantage of the proton wave function we 
find that the interference of the proton spin amplitudes can provide necessary 
phase for Sivers function of down quark. The inclusion of vector diquark has 
also been adopted in Ref. [17], where a version of the axial- vector spectator 
model [18] was employed. Although there are vector diquarks in both mod- 



els, we will explain later that, our model is different from that in Ref. [17]. 
Both models can produce qualitatively similar Sivers functions and asymme- 
tries, but our model has fewer parameters. Relying on Sivers functions of up 
and down quarks we give some prediction for the Sivers asymmetries in vr"*" 
production in SIDIS, and also rough estimation of those in n~ and 7r° pro- 
duction which is currently analyzed in HERMES Collaboration. One feature 
of our model which will be shown is that if we take away the vector diquark 
spectator, the result of our model returns to the result appearing in Ref. [10]. 

This paper is organized as follows. In Sec. 2 we derive the light-cone wave 
function for the two-body Fock state of the proton including vector diquark 
spectator utilizing the Melosh-Wigner rotation. Then we estimate the Sivers 
functions of up and down quarks in the presence of final state interactions in 
Sec. 3. In Sec. 4, we analyze the single-spin Sivers asymmetries of 7r+, vr" and 
TT^ off proton target respectively. Finally we give a brief summary in Sec. 5. 



2 The quark-diquark decomposition of proton in Fock space 



The light-cone formalism [19] has remarkable properties in describing compos- 
ite state such as hadron. It is suitable to describe the relativistic many-body 
problem and there have been many successful applications of the light-cone 
quark model to various physical processes [20]. In light-cone formalism, the 
hadronic wave function is expressed in terms of a series of light-cone wave 
functions in the Fock state basis. For the proton a convenient basis is the 
quark-diquark two-body state, which has been adopted in Ref. [10], where the 
two-body state wave function was derived from the relativistic field theory 
treatment by calculating the interaction vertex in light-cone frame [10]: 



u{k+,k-,k±) u{P+,p-,P± = Oj 



Vk+ VP+ 

u{k~^, k~, k^) and u{P^, P^ , Pj_) are the light-cone spinors of the quark and 
the proton respectively. A similar treatment was employed in Refs. [21,22] 
where the electron-photon two-particle Fock state decomposition was obtained 
from the interaction vertex: 

u{k+,k-,k^) , n(P+,P-,P^) 

-7 ■ e 7^= , (2) 



^A;^ 



in which the light-cone spinors represent the electrons and e'^ represents the 
polarization of the photon. 



The same issue can be also analyzed in the light-cone quark model. In this 
model, the light-cone wave function of a composite system can be obtained 
by transforming the ordinary equal-time (instant-form) wave function in the 
rest frame into that in the light-front dynamics, by taking into account rel- 
ativistic effects such as the Melosh-wigner rotation. The equivalence of two 
approaches has been demonstrated in Ref. [23], where both approaches are 
used to derive the pion light-cone wave function and lead to the same re- 
sult. The Melosh-Wigner rotation is one of the most important ingredients 
of the light-cone formalism, and it relates the light-cone spin state \J,X)f 
to the ordinary instant-form spin state wave functions \J, s)t by the general 
relation [24,25,26] 

|J,A)^ = ^f//,|J,.)T. (3) 



The effects due to the Melosh-Wigner rotation have been calculated for the 
nucleon axial charges [26,27], the magnetic moments [27], the nucleon helicity 
[16] and transversity [28] distributions, the quark orbital angular moments of 
the nucleon [29], and recently for the form factors of nucleons [30]. 

We now derive the light-cone wave function for the quark-diquark Fock state 
of the proton in the light-cone quark model, with both scalar diquark and 
vector diquark. The proton wave function in the SU(6) quark-diquark model 
in the instant form is written as [16,31] 

^t.^(gD) = sin^(^y|g\/)^'* + cosOip s\qS)^'^, (4) 



with 



o 

\qSf'^ = S{ud)u^'\ (5) 

where ff-, -IJ. and t, J, label the spin projection J^=±.^ of the proton and J^ =±| 
of the quark, respectively, VsX(lil2) stands for the qiq2 vector diquark Fock 
state with third spin component s^, S{ud) stands for a ud scalar diquark 
Fock state, and lpd represents the momentum space wave function of the 
quark-diquark state with D denoting the vector or scalar diquarks and 9 being 
a mixing angle that breaks SU(6) symmetry at 6* 7^ 7r/4. Here we choose 
9 = IT /A following Ref. [16]. Applying the transformation Eq. (3) on both 
sides of Eq. (4), we will obtain the spin space wave function of the quark- 
diquark Fock state in light-front frame. The proton state with zero transverse 



momentuin P^ = 0_l in the l.h.s. of Eq. (4) keeps identical under the Melosh- 
Wigner rotation: 

^^'^P+, P~, 0^) = ^^/{P^, P', 0^). (6) 

For the r.h.s of Eq. (4), we first consider the spin part of quark-scalar-diquark 
Fock state,which is determined by the spin of the quark Xq'- 

X\qs) - Xq ■ {' ) 

The Melosh transformation on a quark spin state with four-momentum {k^, k) 
is 



xl{T)=LJ,[{k^ + m)xl{F) - k\l{F)], 

xiiT)=uj,[{k^ + m)xi{F) + k\l{F)], (8) 

where Uq = [2k~^{k^ + m)]~2 , fc^'^ = k^ ± ik5_, and k'^ = k^ + k^ = xM, x is 
the momentum fraction of the quark, m is the mass of the constituent quark, 
and Ai is the invariance mass of the composite state, which is approximately 
taken as the proton mass M in this work. Therefore we get the light-cone spin 
wave function of the quark-scalar-diquark part for the J^ = +| proton: 

X\qS){x,k^) =J2CUi^,^±^Jq)x'^^ (9) 



where the component coefficients Cf^i(x, k_L, J^) have the forms: 



Cfi (x, kx, ^)=ujg{xM + m), 

^+i(^,k±,i) = -^,(ki + ikl). (10) 



With the momentum space wave function for the diquark, we can write down 
the light-cone wave function for quark-scalar-diquark Fock state as follows 



|5g(P+,Px = 0x))^=/ ^''^^^^ [^^(x,k^,T)|xP+,k^,T) 

•^ \/x{l — x)167r'^ 

+^^(x,kx,i)|xP+,k^,i)], (ii; 

where 



^|(x,k^,T) = (M + ^V5, (12) 

4(x,k^,l) = -^±±^^s, (13) 

which is consistent with the hght-cone wave function employed in Ref. [10]. 
For the momentum space wave function of the scalar diquark we adopt 

^ ^/^g^ .... 

'^^ M2-(ki + m2)/x-(ki + A5')/(l-x)' ^ ' 



where Xs is the mass of the scalar diquark. The Ug appearing in Eq. (10) can 
be considered to be incorporated into the momentum space wave function ips- 
It is shown that there is a component tps{x, k_L, |) in the proton wave function 
coming from the Wigner-Melosh rotation effect, which does not exist in the 
non-relativistic constituent quark model. Note that each spin configuration 
satisfies the spin sum rule J^ + L^ = +|, where L^ is the orbital angular 
momentum projection. 

Now consider the quark- vector-diquark two-body Fock state in the right hand 
side of Eq. (4). The spin part of the state composed by up quark and vector 
diquark with J^ = -|-| is written as 

Xvxl - ^X^\l (15) 



where xv stands for the spin of the vector diquark. For the spin-1 vector 
diquark, the Melosh transformations read [32] 

xf{T)=u?y\{k^ + \v)\f{F) - V2{k^ + \v)k^xl{F) + {k^fxy\F)l 
x'viT) =uUV2{k+ + Xv)k^xVin + mv + V)A;^ - k^k^hvi^ 

-V2{k+ + Xv)k^Xv\n}, 
Xv\T) =ujl[{k^)\V{F) + V2{k^ + \v)k^xl{F) + {K + \v)\v\m) 

where Ay is the mass of the vector diquark, ky = ky + ky={l — x)M., ky = 
ky_|_ ± iky_|_, here ky_L is the transverse momentum of the vector diquark 
satisfied kyj_ = — kj_, and wy = [2ky{ky + Xv)]- 

Now substituting Eq. (16) and Eq. (8) into Eq. (15), we obtain the spin space 
wave function of the up-quark- vector-diquark state for the Jp = +\ proton in 
the light-front frame: 

X|y.>(x,k^)= Y. C'^d^,k^,Jv,'J',)xiXxt\ (17) 

jz Jz 2 



where the coefficients read 



^2 
^2 

^2 



(x, 


k^, 


(x, 


k±, 


(x, 


k±, 


(x, 


k±, 


(x, 


k±, 


(x, 


k±, 



+1, t) = wyy[-V2{k^ + Ay)(A;+ + m) - V2(A;+ + Xv)V ^ 
+ 1, i) = Wqwl[^{k+ + Ay)ki - v^(/c+ + \v)\k+ + m)], 
0, t) = w,wl{2[{k^ + Ay)fc+ - ki](A:+ + m) - 2(A;+ + Ay)kl}, 
0, i) = w;,«;^2[-(A:+ + Ay)fc+ + k^ - 2(A;+ + Ay)(A;+ + m)]A;^, 
-1, T) = Wqwl[V2{k^ + Av)(A;+ + m) - V2ki]fc^, 



-l,i) 



WgWy 



[y2A;f(A;+ + Ay + A;+ + m)]. 



:i8) 



Similar to the quark-scalar-diquark two-body Fock state wave function, we 
notice that there are also higher helicity components in the case of vector 
diquark. They also appear in the analysis of the form factor of the nucleon [30]. 
With k^^^ in the wave function, those transverse momentum dependent higher 
hehcity components (+1, |), (0, |), (—1, t) and (—1, [) come from the Melosh- 
Wigner rotation and reflect the intrinsic transverse motion of quarks in the 
proton. Each spin configuration satisfies the spin sum rule J^ + Jy + L^ = +^. 
The component (— 1, |) is the higher twist component which can be ignored in 
contrast with other components. From Eq. (18) we write down the expression 
for the light-cone wave function of up-quark-vector-diquark Fock state with 
F = -\--- 



d^k I dx 



|\/n(P+, P^ = 0^))^ = / , \^\{x, k^, +1, T)|xP+, k^, +1, t) 



-i,l{x, k^, +1, i)|xP+, kx, +1, i) + ^J(x, kx, -1, T)kP+, kx, -1, T) 

-V^^(x,k^,0,t)|xP+,k^,0,t)+^^(x,kx,0,i)|xP+,k^,0,i)], 



(19) 



in which 



i,l{x, k^, +1, t) = -^^^rr^^y^ 



4(x,kx,+l,i)=y2(M 



x(l — x) 



1 — X 



)<^v 



^|^(x, k^, -1, T) = V2^^^^(/^y, 

1 — X 

V'^(x, kx, 0, T) = [2(M + -^) - (M + -) 

L 1 — X X 



V^y, 



7/.^(x,k^,0,i) 



1 -h X k^ + ik^ 
1 — X X 



'-'^v-, 



(20) 



where ^y stands for the momentum space wave function, which has the same 
form of Eq. (14), but A^ replaced by Ay. The above components of the wave 
function are similar to the case of the light-cone wave function of \e'i) Fock 
state in Ref. [21,22], except that there are two components with spin projection 
Jy=0, according to the mass of the vector diquark. We should point out 
that Eq. (20) is the parametrization of the spin coupling of the quark-vector- 
diquark state, in which the three spin projections of the vector diquark have 
been specified explicitly. This is different from the spectator model employed 
in Ref. [17], where the spin coupling is parameterized by a covariant nucleon- 
quark-diquark interaction vertex. 

Similarly we can obtain the Fock state of dowra-quark-vector-diquark which is 
as the same as Eq. (20) only with an extra factor — \/2, as shown in following 
schematic form of the proton light-cone wave function with J^ = -|-i: 

^\ = ^\Su)l + ^\Vu)l - l\Vd)l. (21) 



The Jp = ~^ proton wave function can be evaluated in the same way. 



3 Estimates of Sivers function for up and down quark 



The final state interactions are generated by the gauge link term in the gauge- 
invariance definition of the parton distribution funcitons. In model calcula- 
tions [17,33,34,35] of non- vanishing Sivers function, one gluon-exchange ap- 
proximation for the gauge link exponential was adopted. In the language of fi- 
nal state interactions, the nonzero phase comes from the interference of proton 
spin amplitudes, in which different proton spin states J^ = ±| couple to the 
same final state \F), requiring the orbital angular momentum of the two pro- 
ton wave functions differing by AL^ = 1. In previous section, we have derived 
the light-cone wave function including quark-vector-diquark structure. Like 
the quark-scalar-diquark state in Ref. [10], the quark-vector-diquark structure 
also contributes to the proton spin amplitudes, and further the final state 
interactions. Following the calculation in Ref. [10], we calculate all possible 
proton spin amplitudes, by adopting one gluon-exchange approximation: 



Aiii^ -1 T) = V2'-^±^CiK + i^^.2), 

i — A OTT 

Ait^ +1 i) = -V2{M + -\)C{K + i^^.i), 
A{t^ +1 T) = -V2^iy^^C{K + ^^9vi), 



A{t^ T) = (Af + 2-^ - '^)C{K + i^g.i), 

i — A A OTT 

A(it^ S T) = (M + ^)C(h. + i^g.O 



iTT 



a I ™2 



Aiil^ Si) = -"^^^dK + i^^^2), (22) 

for Jp = +1 two-body states and 



A{\^^ -1 T) = -V2{M + -^)C(/i. + 2^<7.i), 

i — A ovr 

^(^^ +1 i) = -V2^1^CiK + i^9.2), 

1 — A OTT 

1 • 2 

^(^^ -1 i) = V2^^^^^C{K + i^^.2), 
^(4- T) = ^(ri - irl)C{K + i^^.2), 

i — A OTT 

A{1^^ i) = (M + 2-^ - ^)C(/i. + i^^.2), 

1 — A A Svr 

^(4- 5 T) = ^^^C7(/.. + i^^.2), 

A OTT 

A{\^^ Si) = {M+ '^)C{K + i^^.i), (23) 

A Svr 

for Jp = —\ two-body states. Where A and r^ are the Bjorken variable and 
the transverse momentum of the struck quark respectively, ei and 62 are the 
charges of the struck quark and spectator diquark, and 



C = -^eiP+VA2A(l - A), (24) 

^'^' ^ rx + A(l-A)(-M2 + f + ^)' ^^^^ 

in which g is the proton-quark-diquark coupling constant. In literature there 
are two choice for g, one is to take g as a constant [10,33], another as a form 
factor [17,35] relying on the momentum of the outgoing quark, to converge 
the transverse momentum integration of the distribution functions. As argued 
in Ref. [33], the inclusion of the form factor would add another complication 
on the evolution of the asymmetry, therefore we treat g as a constant in our 
model, which is equal for the scalar and the vector diquark. The last two 
equations in Eq. (22) and Eq. (23) are as the same as Eq. (5)-Eq. (8) in 
Ref. [10], while the other ten equations are contributed from vector- diquark 



structure. With above amplitudes we calculate the asymmetries in quark level 
for up and down quarks through [11] 



_ im{r:^i') 



qyu/d 



Jl^u/d 



(26) 



Here T^^ and 7^^^ are the interference terms coming from the final state in- 
teractions with the forms 



Vnt = A{t^ S T)*^(4- 5 T) + -A{t^ -1 T)M(4- -1 T) 



^(^^Ot)*^(^^OT), 



%nt = gAt^ -1 T)M(^^ -1 T) + -A{t^ r)*^(^- T), 

respectively, and TW"/*^ are denoted as: 



(27) 



A^": 



1 



A{t^ S T) 
1 



' ^ A{t^ 5 i) ' + ^ 



A{t^ +1 i) 



2 1 



Aiii^ -1 i) 



A{t^ +1 T) 
Ait^ t) 



^ ^(^^Oi)' 



A^' 



9 



2 1 



2 1 



^(^^+1T) +g^(^^+li) +-A{t-^-ll) 

^(^^OT)f + ^|^(^^Oi) 



(2J 



Substituting Eq. (22) and Eq. (23) into the equations above we obtain 



r.u -_lfM + -)i^i^ ^ In^^^ 



r2 



(M 



Xi 



1- A'l- A 9 



"^ +l(M + 2 ^'' 



m. 



1-A A- 



1-A 



6162 



In 



A.(r^ 



1 + A) AJ Svr A.(ri)ri A,(0) 



%: 



6162 



int ' 



In- 



A„(r2 



47r A„(ri)ri A,(0) 



Xi 



9^ 1-A'l-A 



4'^' + 2r^A A'(1 + A)A 



m 1 — A r-*; 
a"/ 



(29) 



(30) 



10 



M"" 



ra. 



(^ + ^r + ^. 



hi + 



9 



(M + 



A 



V n2 



1-A' (1-A)2 






m 



A2(i_A)2 (1 + A)2 2A2 ^2 1-A 2A- 



M'^=- (M + 



A 



V n2 



+ 



+ 



(1 - A)2 ri 



9L^"" ■ 1-A' (1-A)2 A2(l-A)2 (1 + A)2 2A2 
Ay m 



"2 1-A 2A' 



h. 



(31) 



(32) 



in which 



K/v{vl) = vl + A{l-A){-M^ + 



m 



A 



5/y 



(33) 



From Eq. (26) and the relation between the asymmetries Vy, f^{A,r 
/i(A,ri) [11,33] 



and 



pa 



y -^/iT(A,rl)/A"(A,rl), 



(34) 



we obtain the Sivers functions of up and down quarks by 



-pa 



/ii,'^(A,ri) = --f/r(A,ri). 



(35) 



Here a is the flavor index including u and d. The difference between Eq. (34) 
and Eq. (32) in Ref. [33] is that we assume that the equation stands for both 
up and down quarks. From Eq. (34) we can see that if we exclude vector 
diquark, the result of our model will return to the results in Ref. [10] and 
Ref. [33]; when vector diquark is included, the Sivers function for down quark 
is evaluated. There is another advantage of Eq. (26): the coupling constant g 
will appear in both the numerator and denominator of Eq. (34) and cancel, 
thus the Sivers functions in our model is free from g and normalization factors. 
To avoid ambiguity, we use below the same notation used in Ref. [10] with A 
denoting the Bjorken variable of parton and r^ for the transverse momentum 
of quark. In realistic application the moment of the Sivers function is more 
useful. The ri-moment of the Sivers function is defined as [36] 



A-«(A)^/.V(^)Ai(A, 



(36) 
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Fig. 1. The r^-moments of Sivers functions of u and d quarks in SU(6) light-cone 
quark-diquark model. 

In numerical estimates of this moment we adopt Gaussian distribution of the 
transverse momentum for the unpolarized distribution function 



/i(A,r^ 



/i(A) 



Exp 



{r±y 



(37) 



in order to guarantee convergence of the r^ integral in Eq. (36). Here (r^) is 
the average transverse momentum of the quark. For the transverse momentum 
integrated unpolarized distribution /i(A) we adopt CTEQ6M parametriza- 
tion [38] for the valence quarks. Also we choose the parameters in the expres- 
sion as: \s = 0.6 GeV, Ay = 0.9 GeV as estimated to explain the A^-A mass 
difference, m = 0.36 GeV, M = 0.94 GeV, (r^) = 0.4 GeV. To generalize the 
analysis to the corresponding calculation in QCD, one needs to extrapolate 
\^^\ — > Cpds with Cp = |. We choose as = 0.2 to correspond the kine- 
matics of HERMES experiment. With above parameters we obtain the Sivers 
functions of up and down quarks explicitly. The r^-moments of the Sivers 
functions of up and down quarks are figured by Fig. 1 which shows that in our 
model the Sivers function of up quark is positive, while that of down quark 
tends to be negative, and the former is 6~7 times larger than the later in 
scale. This agrees with the r^-nioments of the Sivers functions that appeared 
in Ref. [17]. And the size of the moments is similar to that of the bag-model 
result presented in Ref. [34]. 



4 Single-spin Sivers asymmetries of pion electroproduction in SIDIS 



Sivers effect can be observed in meson electroproduction in semi-inclusive DIS 
with unpolarized electron beam off transversely polarized nucleon target which 
is measured in HERMES. With the Sivers functions of up and down quarks 
obtained above, we calculate single-spin asymmetries in SIDIS contributed 



12 



from Sivers effect by selecting proper weighting factor J-^sin(0 — (ps)- is 
the azimuthal angle between the lepton scattering plane and the transverse 
momentum of outgoing hadron, (ps is the azimuthal angle of the spin of proton 
target. The asymmetries are expressed by the convolution of the r^^-moment 
of Sivers function and the usual unpolarized fragmentation function Di[z), 
divided by the unpolarized cross section [37] 

l^-U(,-«\ (A..)^?#C^1^. (38) 



M ^'^ ^"/vT ' E„e^/f(A)Z);(c) 



where the subscript UT represents unpolarized beam on transversely polarized 
target, Ph± is the transverse momentum of the pion. For Di{z), we consider 
both the favored fragmentation function 

D{z) = Df{z) = Dl{z) (39) 



and unfavored fragmentation function [7] 

D{z) = Df{z) = D:-{z) (40) 

for TT^ fragmentation. Also we assume 

D^'^iz) = l[Df{z) + D:-{z)] = l[D{z)+D{z)] (41) 

for vr" fragmentation. For the explicit analytically forms of D{z) and D{z) we 
adopt [39] 



D{z) = 0.689^-^-°=^9(l - zy-^^\ 

D{z) = 0.217^"^-^°^(1 - z)2-°^^. (42) 



For 7r~ production, the favored fragmentation process is d -^ it~ . Since we 
have obtained the Sivers function of down quark, we can give some predic- 
tion of the azimuthal asymmetries for the vr" production. We have ignored 
the sea quark contribution, and assume that the contribution comes mainly 
from the valence part, since we can not get information about sea quark from 
the SU(6) quark-diquark model. For comparison with future HERMES exper- 
imental data we also calculate the asymmetries averaged by the kinematics 
of HERMES: 0.023 < x < 0.4 and 0.2 < ^ < 0.7. The single-spin Sivers 
asymmetries of tt^, tx~ and 7r° plotted vs x and z are presented in Fig. 2. The 
figure shows that the x-dependent asymmetries of 7r+ production are larger 
than that of vr" production. The 2-dependent asymmetries of vr"*" and 7r° are 
similar, and different from that of vr" in shape and size. 
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Fig. 2. Single-spin Sivers asymmetries measured in SIDIS off proton target. The 
asymmetries liave been weiglited by -^j- The first, second and third rows corre- 
spond to the asymmetries of vr"^, vr" and vr*' production respectively. The left column 
corresponds to the x-dependent asymmetries, the solid and dotted curves corre- 
spond to the asymmetries averaged and not averaged by the HERMES kinematics 
0.2 < z < 0.7. The right column corresponds to the z-dependent asymmetries, the 
solid and dotted curves correspond to the asymmetries averaged and not averaged 
by the HERMES kinematics 0.023 <x < 0.4. 

5 Conclusion 



Final state interactions provide a reliable explanation for Sivers single-spin 
asymmetries in SIDIS. To get phenomenological estimation we derive the light- 
cone wave function of the proton based on the SU(6) quark-spectator-diquark 
model. We obtain the light-cone wave function containing both scalar diquark 
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and vector diquark. The transverse momentum dependence, which is responsi- 
ble for the transverse momentum distribution function and single-spin asym- 
metries, is introduced to the wave function by employing the Wigner-Melosh 
rotation. The Wigner-Melosh rotation also brings the higher helicity com- 
ponents which contribute to the final state interactions. Utilizing the wave 
function we calculate the Sivers functions of up and down quarks in one-gluon 
exchange. We find that in our model the Sivers function of up quark is pos- 
itive, while that of down quark tends to be negative, and that the former is 
much larger than the later in scale. We give analysis of —^ weighted Sivers 

asymmetries A^^^ of vr^, vr^ and 7r° production off transverse proton 

targets plotted vs x and z from our model. The numerical result shows that 
the x-dependent asymmetries of ir'^ production are a little larger than that of 
7r~ production. The 2;-dependent asymmetries of vr"^ and vr*^ are similar, and 
different from that of vr" in shape and size. 

In contrast to the result in Ref. [17], our Sivers functions and asymmetries 
are qualitatively similar but quantitatively different. We conclude that the 
differences come from two aspects. On the one hand, we use different quark- 
vector-diquark structure from that of Ref. [17]; on the other hand, we treat 
proton-quark-diquark coupling g a.s a. constant, different from that of Ref. [17]. 
As a result, the Sivers functions in our model are free from g and normalization 
parameters, but rely on (r^). 

Similar to other model calculations, we use some assumptions and approxima- 
tions, such as one Abelian gluon approximation and valence quark dominance, 
and do not consider the evolution of the Sivers function, thus at current level, 
our Sivers functions and asymmetries are rough estimates. 
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